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New corrections to the axial vector currents coming from the sub- 
leading terms in the 1/N C expansion are shown to be non-vanishing and 
substantially improving the phenomenological predictions for axial quan- 
tities. These corrections are shown to come from two distinctive sources: 
1) anomalous part of the Euclidean effective action related to the Wess- 
Zumino term of the SU(3) Skyrme model and 2) real, non-anomalous 
part which in this order of 1/N C has no counterpart within any local ef- 
fective meson theory. The appearance of the type 2) terms is due to some 
explicit time- ordering of the collective operators entering the formulae for 
the axial constants. The question of regularization of these quantities is 
discussed. The analytic symmetry breaking terms in the strange quark 
mass play a minor role for and g^ . They are however important for 

Finally the numerical values for the <7a's are g^ = 0.37, g^ = 1.38 

(s) 

and g A = 0.31 reproducing reasonably well the recent data from lepton 
scattering. 
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I. INTRODUCTION 



It is a long lasting problem to determine static properties of hadrons from the gen- 
eral theory of the strong interaction, Quantum Chromo Dynamics (QCD) [|J. There- 
fore there were some attempts in the past to derive an effective theory for the strong 
interactions in some low energy approximation . However none of the derivations 
could exactly claim the range of validity of the approximations. Although the resulting 
theory, which coincides with the formerly invented Nambu-Jona-Lasinio (NJL) model 
J7U9], does not confine it shares the maybe most important features of low-energy QCD 
relevant for ground states. These are: chiral symmetry and its spontaneous breaking. 
In this case the lagrangian itself is chirally invariant but the vacuum state breaks the 
symmetry. The symmetry breaking is driven by a non- vanishing quark condensate, 
leading to the constituent masses of the formerly massless quarks. As a consequence 
the Goldstone bosons emerge, namely the pion, kaon and eta. These almost massless 
excitations of qq pairs are expected to play a dominant role in the QCD vacuum not 
only in the long range limit but also at small distances P1,|T0 . 



In the presently investigated NJL model the nucleonic scenario is realized by these 
Goldstone bosons and explicit quark degrees of freedom. The gluon degrees of freedom 
are already implicitly contained in the lowest order because the action is obtained, 
at least formally, after path integration over the gluon fields. The quarks are then 
bound in a selfconsistent potential based on a non-trivial chiral field configuration in 
the Hartree approximation [|Tl|-|T3| which is leading in 1/N C expansion. Though this 
picture resembles much an effective meson theory like the Skyrme model [H ], one 



should note that here the soliton is non-topological, i.e. practically speaking based 
on the non-perturbative dynamics. On the contrary the Skyrme soliton is topological, 
i.e. the baryon number totally hinges on the topological winding of the chiral field, 
while in the NJL model the baryon number is carried by 3 (or N c ) valence quarks. On 



the basis of the analysis of the proton polarization it was argued in Ref. |15| that the 
valence quark degrees of freedom, missing in the Skyrme model, do account for some 
important physics. 

There is however another, more profound difference between the present model 
and the effective meson theories. The effective meson theories are based upon some 
lagrangian density being local in time and space, whereas the NJL type models are 
formulated in terms of the path integral over the fermions, which can be regarded as 
a time- ordered product of the field operators. It has been for a long time overlooked 
that this time-ordering may bring up some new contributions for various baryonic 
observables. In the recent papers [p^ -[T8|| these contributions (hitherto referred to as 
time- ordered) have been calculated for axial decay constants and for magnetic moments. 

The example of the axial constants is perhaps the most persuasive. It is well 
known that in the nonrelativistic quark model g& = (N c + 2)/3. This means that 
there are important 1/N C corrections to g&, which for N c = 3 amount to 60% of the 
leading result. In the effective meson theories the leading term for gx scales also 
as N c , however the next-to- leading correction comes only at the 1/N C level in the 
SU(2) version of the model. This is due to the fact that effective meson lagrangians 
are even in field derivatives. In the cranking approximation for the rotating soliton 
each time derivative counts as 1/N C . The only possible source for the contributions 
linear in the time derivative is the Wess-Zumino term which vanishes identically in the 
SU(2) case. However in the SU(3) case both in the NJL and in the Skyrme model 
such corrections appear. In the NJL model in the leading order of LWLA (Long 
Wave Length Approximation) or gradient expansion they are equal to the Witten's 
anomalous current 0] of a local mesonic theory. 

It was always believed that the effective quark theories are equivalent to the effective 
meson theories in a sense that they correspond to some calculable local lagrangian 
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density. Strictly speaking such a lagrangian density can be derived from the NJL like 
model only in the limit of the large soliton size. Although this statement remains true 
it does not mean that the matrix elements of some operators, take as example the 
axial currents, are the same irrespectively if they are calculated straight away from the 
fermion path integral or from the equivalent local lagrangian density. This is due to the 
two facts: 1) upon the semiclassical quantization of the soliton the cranking velocities 
are promoted to the collective operators which do not commute with the rotation 
matrix itself and 2) the path integral is time- ordered, i.e. it dictates unambiguously in 
which order the cranking velocity and the rotation matrix appear in the expressions 
for the matrix elements of the axial currents. If these non-local properties of the path 
integral are properly taken into account then one gets the desired 1/N C corrections. 
These corrections are not small and improve the phenomenological predictions of the 
NJL model. On the contrary in the local limit of the present effective quark theory 
they are identically zero. 



In the previous paper [17j we have calculated the three SU(3) axial decay constants 
9a\ 9a^ an d 9a in the chiral limit within the semibosonized NJL model. This model 



reproduces the hyperon spectra []20| and also the isospin splittings within baryon mul- 
tiplets ||21|| . Earlier the properties of the axial currents have been investigated in the 
NJL model only for the case of SU(2) P2]JT3[ ]. Then the 1/N C corrections for SU(2) 
have been roughly estimated in |T(| , neglecting regular izat ion and sea contribution and 
to full extent in fl8|]. Beyond this there are only calculations of the axial coupling con- 



stants within the pseudoscalar SU(3) Skyrme model |23[ and the pseudoscalar vector 
meson SU(3) Skyrme model |24j| . Actually hyperon spectra indicate [^(J that the scalar 
and pseudoscalar NJL model gives a more refined structure of the collective hamilton 
operator, than the pseudoscalar Skyrme model. A comparable structure can be ob- 
tained in the Skyrme model only by introducing explicit vector mesons and modeling 
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the anomalous and symmetry breaking part of the effective action. This introduces a 
large number of parameters, whereas in the NJL model the parameters Q can be fixed 
completely by requiring proper mesonic masses and decay constants. The symmetry 
breaking pattern and the anomalous effective action is then uniquely determined. 
In contrast to our previous paper |T7| now we implement explicit time- ordering 



within the framework of the effective Euclidean action (EEA). This treatment allows 
us to make a clear distinction between the terms which emerge from the real or form the 
imaginary (anomalous) part of the EEA. Actually it turns out that the new explicitly 
time- ordered terms emerge from the real part of the EEA and therefore have to be 
regularized. In the present approach the regularization prescription is unique and the 
regularization function is derived in a well defined manner. 

In the present paper we furthermore extend our previous calculations and calculate 
the m s corrections to the axial coupling constants. Besides the theoretical interest 
due to the new explicitly time- ordered corrections, the axial constants are of utmost 
phenomenological importance as far as the comparison with the recent measurements 
of the polarized proton and neutron structure functions is concerned. So the main 
phenomenological concern of this paper will be a comparison of the model predictions 
with the experimental data for these quantities. Especially we will concentrate on 
the role of the explicitly time- ordered corrections and furthermore on the so called 
anomalous quantities, which are dominated by the valence contributions [ERE^jnJ]. 
The latter manifest another conceptual difference between the Skyrme model and the 
present NJL model. 

The organization of the paper is as follows. In Sec. II we review the basic features 



1 The free parameters are: M - constituent quark mass, and cut-off parameters and to some 
extent m s - strange quark mass 
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of the NJL model with special emphasis on the solitonic description. In Sec. Ill 
we describe the quantization procedure and summarize the results on the hyperon 
splittings. We use the mass splittings to fix the parameters of the model. In Sec. IV 
we derive expressions for axial currents in the chiral limit. Special emphasis is put on 
the new contributions from the explicit time-ordering and their regularization. Then 
in Sect. V we discuss mass corrections to the axial currents. Our numerical results are 
presented in Sect. VI. Section VII contains a brief comparison with the results of the 
Skyrme model. We present our conclusions in Sect. VIII. 

In the Appendices we present useful formulas for the semiclassical quantization 
(App. A), a derivation of the regularization functions (App. B), a gradient expansion 
for normal and time-ordered quantities (App. C). 

II. THE SU(3) NAMBU-JONA-LASINIO MODEL - SOLITONS 

The quark Nambu-Jona-Lasinio model |7],|| can be written in the four-fermion 
formulation as: 

£njl = q(x)(ip - m)q(x) - 2G[(q\ a q) 2 + (qij 5 \ a q) 2 }. (1) 



Here the summation over the A a matrices is implicit (with Ao = y 2/3) an d Tn is the 
bare quark mass matrix. In the chiral limit the Lagrangian has the desired SU(3)r® 
SU(3)l symmetry in addition to the U(l)v® U(i)a, where the U(1)a is the symmetry 
which is not shared by QCD. In principle one could introduce the 't Hooft term into 
the Lagrangian, which breaks the U(1)a explicitly. It could then serve as a source for 
the r)' mass, which otherwise would be a Goldstone boson. This was recently done 



by Kato et al. |26j and the resulting profile for the SU(2) soliton was very similar to 
the solutions that are restricted to the chiral circle (i.e. non-linear case) and that are 
used here. So we conclude that the effects of the 't Hooft determinant on the solitonic 
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observables of the present calculations are rather small. However it was proven in Ref. 



[ 26 1 that the 't Hooft term leads to a stabilization of the linear version of the model, 
which furthermore supports the use of the non-linear Ansatz for the chiral fields. 



Performing the bosonization procedure of Eguchi [27| one arrives immediately at 
the new classical Lagrangian 

£e = q(x)(-ip + g (a a X a + i l5 n a \ a ) + m)q(x) + )y 2 (a a a a + n a n a ) , (2) 

where the original coupling G is given by G = g 2 /fi 2 - For convenience Eq.(||]) can be 



rewritten in the polar decomposition of the chiral fields |28] : 

1 u 2 

£ E = ?(z) + PkM + PlM + + m) q{x) + — ^-Tr^ (lM f ) , (3) 

where Ai = g(cr a X a + i r y^TT a X a ) = £lM£r is given in terms of the unitary matrices £l, £r 
and a hermitian matrix M. The -Pr/l = l/2(l±7s) are the right/left helicity projection 
operators. Furthermore the unitary gauge £r = = ^ can be chosen to eliminate the 
redundant degrees of freedom. As an approximation the scalar degrees of freedom in 
M are frozen, such that M is just the constituent quark mass matrix and £ is related 
to the chiral field by U = £ 2 . The latter can be parametrized as exp {m a X a / f^). This 
parametrization corresponds in SU(2) to the chiral circle condition a%\ +if 2 = const, 
where 07 2 ) is the SU(2) isoscalar a field. Then the last term in eq. (|[) is just a constant 
and will be omitted in the following. 

Now we will shortly summarize how to calculate quantum corrections to Ce which 
are due to the fluctuations of the fermion fields. The fermion determinant obtained by 
the integration over the quark fluctuations corresponds to the leading contribution in 
1/N C expansion. The next to leading terms, the so called boson determinant, have been 



2 We will always denote by E Euclidean and by M Minkowski quantities if it is not obvious 
from the context. 
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shown to have minor influence on the mesonic observables | 29fl . Hence the effective 
Euclidean action (EEA) can be written as: 

S eS [U(x)} = -Sp log (-ip + P R M(x) + P L M^(x) + m) , (4) 

where Sp is the functional trace in color, flavor and momentum space. 

From eq. (^) the parameters of the model can be fixed by requiring experimental 
values for the pion decay constant f n = 93 MeV, the pion and the kaon mass, m n = 



139 MeV and mx = 496 MeV, (see Ref. [2(J for details). As a result the constituent 



quark mass is the only free parameter of the model, which can e.g. be used to fix 
baryonic properties pjjffiH . 

Solitonic solutions of eq. (^) can be found by making a time-independent hedgehog 
Ansatz for the chiral field U and writing the SU(2) action in the chiral limit in terms 
of a single particle intrinsic hamiltonian: 

H = -i 1 ±{-i li d i + MU{x)). (5) 

We will specify the SU(3) extension of the SU(2) Ansatz in the next Section. In the 
proper time regularization ||30|| , the effective action becomes in this case: 

S eff = Sp J^ u )exp(-u(-d 2 T + H 2 )), (6) 



where <p(u) = c8(u — 1/Af) + (1 — c)8(u — 1/A|) is the regularization function of Ref . [21J, 
which reproduces common values for the current quark masses and quark condensates 
in the vacuum. The classical equations of motion can be solved selfconsistently for the 
chiral field U, resulting in a localized soliton with unit winding number. The energy 
spectrum of the Hamilton operator H (eq. (§)) for the baryon number one sector 
contains a discrete valence level inside a mass gap of the size 2M []12 , |ll ,31|. Then the 



classical energy of the soliton can be written as [^T 



M cl = A c £ val + N c E sea _, (7) 



where E va \ is the energy of the valence level and E seSL resembles the polarization the 
Dirac sea as a sum over the whole spectrum of the Hamilton operator H. 



III. QUANTIZATION OF ZERO MODES AND MASS SPLITTINGS 

The purpose of this section is to apply the semiclassical quantization method to 
the solitons |T!^]TT| . |3"2f| of Sect. [IT], which result from the classical and time- independent 



equations of motion. The idea hereby is the following. In order to quantize the system 
one can perform a time-dependent transformation [Q, which can be in the direction 
of the symmetry or orthogonal to it. If the symmetry is at least an approximate 
symmetry then excitations in this direction should be the dominant contribution to 
the low lying resonances of the model. In order to check this numerically in the present 
model, we consider in addition to the usual expansion of the EEA up to the second 
order in the rotational velocity consistently the quadratic corrections from the strange 
symmetry breaking terms. 

Therefore following the treatment of Ref. |I3| we quantize the soliton by introducing 
time-dependent SU(3) rotations and impose canonical quantization conditions for the 
collective coordinates of the rotation matrix. This will allow for the definition of 
generators of SU(3) and the corresponding baryon states. 

First we make use of the trivial embedding of Witten j31| of the SU(2) chiral field 
Uq(x) = ((7(2) + ijsTTt)/ fn into the isospin subgroup of SU(3) according to 

U 0' 
1 

The soliton solutions of SU(2) are also solutions for SU(3). In the quantization proce- 
dure the embedding (§) generates, as will be seen later explicitly, the correct quantum 
numbers for baryons ||34|| . 

Next one introduces a time-dependent rotation: U(x,t) = A(t) U(x) A(ty. This 
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U{x) 



rotation can be undone by rotating the quark fields: q = A{tyq and q = qA(t). Then 
we obtain: 

AU = ift E = \K^l (9) 

and the relation between Euclidean and Minkowski velocities holds: iQ E = and 

Expanding Sl^ = — Sp \og(d T + H+iQ E — ij^A^mA) of Ref. [20[] up to the quadratic 



order in Q (in Minkowski metric and in the chiral limit) one gets: 

l = -m c1 + hi a i ab n b - ^=n 8 (io) 

where tensor of inertia I ab = diag(/i, I\, Ii, I2, h, h, h, 0) can be found in Ref. |20| . 



The original path-integral / T>U(x, t) will be in the following approximated by the 
integral over the rotation matrices A(t) only, neglecting translations and other fluctua- 
tions |3l| . This is known as the quantization of the rotational zero modes |33| . Instead 



of functionally integrating over all A(t), one makes usually use of the fundamental 



J VA{t) exp {-J dtl%) =< A(T/2) \ exp(-TH (0) ) \ A(-T/2) >, (11) 



relation between the path- integral and the hamiltonian operator formulation |3q] : 

r T/2 
l-T/2 

where is a collective rotational hamiltonian corresponding to the Lagrangian of 
eq. fllOD . In this way the path- integral can be evaluated in terms of the eigenstates of 
the collective hamiltonian. When we consider expectation values of currents similar 
relations hold. There however one has to pay attention to the fact that the path- 
integral is time-ordered in a natural way [ p5| , p6[ , whereas in the operator formalism 
time- ordering has to be introduced explicitly. The latter fact will be of importance 
when we consider expectation values of currents in the next section and is extensively 
described in the case of SU(2) in Ref. Before we calculate axial properties we 



shall concentrate on mass splittings. This allows to judge the perturbation expansion 
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in m s and to fix the remaining parameter of the model - the constituent quark mass 
M. 

To calculate the mass splittings one has to expand the effective action in powers of 
the current quark mass m = [Aq Ao — /is As — ^3 A3 with 

111 
fi = —7=(m u + m d + m s ), |i 8 = -=(2m s -m u -m d ), fj, 3 = -Am (12) 
VD V12 ^ 

where Am = — m u . 

Here an important remark is in order. There are apparently two small parameters 
in the present approach: 1/N C and m s . Unfortunately from the explicit calculations 
one cannot deduce what actually sets the scale for the m s corrections. In any case 
1/N C and m s can be treated as being of the same order. Therefore to be consistent we 
expand the effective action up to terms of the order of m s , m 2 , m s Q and Q 2 (expansion 
in Q corresponds to expansion in 1/N C ): 

L m = -am s + am s D { ^, (13) 
L m n = --y=m s DfjK ab n b , (14) 
L m2 = 2 -ml{N Q {l - D®f + 3N ab D^D®), (15) 

where the constant a is related to the sigma term E = 3/2 (m u + mj) a and 
D% = 1/2 Tr(A^\ a A\ b ). Therefore we can define in this order L tot = Lq + 
L m + L m Q + L m 2. The mass spectrum obtained with the help of L + L m + L m n 



was discussed in Refs. f20| , |21|1 ; there one can also find explicit formulae for K ab = 
diag(-R'i, Ki, K\, K2, K2, K2, K2, 0). Let us here only remind that the anomalous mo- 
ments of inertia are nearly entirely given by the valence part, whereas the con- 
tribution of the valence level to Jj amounts to approximately 60%. The new fea- 
ture of the present calculation is the presence of the moments of inertia N ab = 
diag(iVi, N x , N u N 2 , N 2 , N 2 , N 2 , N /3) in L m 2 defined as: 
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N ab = ~r < 171 I Vft \ nxn \ \ b j Q | m > K p (E n , E m ), (16) 



4 

n,m 



where TZp(E n , E m ) is given by 

1 r du 



n (E n ,E m) = -J^(u 



E n 6 71 E m 6 



(17) 



E — E 

J -'m 

which differs from the regularization function for the usual moment of inertia 
lZj(E n , E m ) |20| because of the different hermiticity behavior of the mass term and 
the Coriolis term (f2) in S^jj . 

The values of iVo,i,2 together with the values of Ji^ and K\ 2 for different constituent 
masses are listed in Tab. |. 

The Lagrangian of eq. (|TUp and eqs. (|TB|-|T4]) reminds the Skyrmion Lagrangian with 
vector mesons (c.f. Ref. P|j|). The quantization proceeds as in the Skyrme model; one 
defines the quantities (see App. A for details): 

N c 

Ja = —Ra = lab^b — ^iDibKb a — (5 ffl g ^^ (18) 

(i = 3 and 8, a, b = 1---8) which are promoted to the spin operators J a = —R a - 
The flavor operators read: T a = —D a i,Jb. Note that despite the fact that J a fulfil 
the SU(3) algebra, only Ji,2,3 have the meaning of the symmetry generators. That is 
due to the structure of the SU(3) hedgehog Ansatz and is reflected in the fact that 
Js = —N c /\/l2 generates a constraint. Therefore the wave function of the baryon 
state B = Y,T,T 3 , J, J 3 belonging to the SU(3) representation 1Z reads (see App. A): 



TZ,B> = VdimJl(Y,I,I 3 | D {n) {A) \ -Y',J,-J 3 ) , (19) 



where the right hypercharge Y' is in fact constrained to be —1. The lowest SU(3) 
representations which contain states with Y — 1 are: 71 =8 and 1Z =10. The quantized 
collective hamiltonian H tot from 

Q J tot 
a ^ 

13 



reads: 



#(°) = M cl + H t 



SU(2) + #SU(3) 



(21) 



Hs\J(2) 



1 



C 2 (SU(2), H sv{3) 



1 



C 2 (SU(3)) - C 2 (SU(2)) 




2/i 



2/ 2 
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Here C 2 denote the Casimir operators of the spin SU(2) and flavor SU(3). M c \ is the 
classical soliton mass. It has been calculated by many authors and its value turns 
out to be relatively large: M c \ rs 1.2 GeV. This is a common problem for all chiral 
models. There are however some negative corrections to it, like Casimir energy or 
rotational band corrections which might bring M c \ to the right value. In this paper, 
instead on insisting on the calculation of the absolute masses, we will concentrate on 
the mass splittings. These are determined in the present model by analytic strange 
mass contributions in 0(m s ) and 0(m 2 s ). Whereas linear terms are given by eqs. 
dl3|Jl4]) , the the various contributions from quadratic m s corrections can be classified 
into the following three cases: 

• quadratic terms from the expansion of the EEA, corresponding to the Lagrangian 
in eq. (p75|) . This will be referred to as kinematical correction. 

• quadratic terms from replacing the rotational velocities in eq. ( |IBD by the gener- 
ators J a by using eq. (j2"Up. This will be referred to as dynamical correction. 

• quadratic corrections from sandwiching the linear part L m ,L m Q, which can mix 
various representations of SU(3), with linear mass corrections from the collective 
wave-function, as will be discussed in the following. This will be referred to as 
wave-function correction. 

Then the hamiltonian from eq. ( ^0[) up to terms quadratic in m s reads: 
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ffffi = I |r 2 X 2 (l - £&) + (ntfi - r 2 K 2 ) £ J m 2 , 



(22) 



H (2) _ 

-"dyn - 



(iV + 3A 2 ) - 2N D 88 + (iV - 3N 2 )D 2 88 + 3(7^ - N 2 ) ]T D 2 8A m 2 s , 



A=l 



where r\ = Ifj/ij. According to items above we have split the 0(m 2 ) hamiltonian into 
the kinematical part i?^, and the dynamical part .£Qy n - 

The Hamiltonian if ^ mixes states of different SU(3) representations, therefore the 
wave function is no longer a pure octet or decuplet but rather a mixture: 

| B > =\S,B > 110,5 > +cf 7 (27,5 >, 

| B' > = | 10,5' > +cf 7 ' | 27, 5' > +cg | 35, 5' > (23) 

where 5=N,A,E,S and B'=A, £*, S*, fi. The coefficients c£ depend linearly on m s , 
therefore with this accuracy there is no need to change the normalization of the wave 
function. In the following we will need their explicit form only for the octet-like states: 



b VE. 
cjq = —(cr-n; 
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I 2 m s , c 27 = — (3cr + n -4r 2 ) 



3 
2 

V6 



(24) 



in the basis [N, A, E,S]. In Fig.l we plot cj| in dependence on the constituent mass 
M. The corresponding 0(m 2 ) contribution to the energy reads: 

,(2) 



E, 



wf 



-{ — (l + Y -X + -Y 2 
\60 V 2 



L fl§ + h - X - -Y 2 ] -(3a + n - 4r 2 ) 2 l J 2 m 2 
250 V 2 2 4 / 9 V ; J 



(25) 



for the octet and for the decuplet: 

5S = -{l(l + ^ + V)l(3.-5, 1 + 2, 2) 2 

l"^-^ 2 ) (* + r 1 -2r 2 ) 2 }l 2 m 2 . 



336 



(26) 
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Here X = 1 — 1(1 + 1) + 1/4 Y 2 is the usual combination entering Gell-Mann-Okubo 
mass relations. 

With the help of the matrix elements of the D functions and spin operators dis- 
cussed in the App. A one arrives at the following result for the hyperon splittings: 

AM (8) =A _L Y -^=X-GY 2 , 
2 v/5 

AM (10) = B - -^-= Y - H Y 2 . (27) 
2y 2 

Constants A and B do not contribute to the splittings within the multiplets, however 
they shift the mass centers and contribute to the 10-8 mass difference. Constants G 
and H not present in the first order Gell-Mann-Okubo mass formula are of the order 
of m 2 . Experimentally one gets: 

F = E - N = 379 MeV, 
D = JL_(E-A) = 86 MeV, 

G = -(3A + S) — -(N + S) = 6.75 MeV (28) 
4 2 

for the octet. For the decuplet the three operators: 1, Y and Y 2 do not form a complete 
basis and therefore there are two independent relations which determine constants C 
and H with some small uncertainty: 

C = V2(E* - A) = -^(Q - 2A + £*) = 422.5 ± 3.5 MeV, 

H = i(2£* - S* - A) = i(3£* -2A-Q) = 2.83 ± 0.33 MeV. (29) 
2 6 

In Tab. |II] we list the coefficients A . . . H for a typical value of m s = 180 MeV as 
functions of the constituent mass M. It can be seen that in order to reproduce the 



experimental numbers of eqs.(|28|,|2"9"|) one has to take the constituent mass of the order 



of 400 MeV. Then all constants A. . . H are roughly reproduced. The constant G and 



H being of the order 0(m 2 ) are small. For reasonable strange quark masses 0(m. 
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corrections to A, B, C and D are of the order of 20% of the leading 0(m s ) terms with 
the exception of F for which 0{m%) corrections are almost zero. 

In order to make phenomenological statements we adopt the following procedure: 
first for given M we find the optimal m s which reproduces 10-8 splitting. To this end 
we define the mean octet and decuplet values: M 8 = 1/2 (A + E) = 1155 MeV and 
Mio = £* = 1385 MeV. Then Ai _b = M 10 - M 8 = 230 MeV is given by: 

Aio- 8 = ry- + B - A. (30) 
All 

Since A — B = const, x one can numerically solve eq.(^) for m s . The result is 
plotted in Fig.2. 

In Fig. 3 we show the m s dependence of the deviations theory - experiment for each 
hyperon. On should remember that for each m s the optimal constituent quark mass 
M was used, so that A 10 _ 8 was automatically reproduced for each m s . The smallest 
deviations ±7 MeV for all splittings correspond to m s pa 185 MeV, i.e. M pa 425 MeV. 

It is interesting to examine to what extent the new corrections calculated in this 
paper are important. The Yabu-Ando method of diagonalizing the hamiltonian of 
0(m s ) exactly is widely spread in the literature |37| , |24]^0| , |38[| and it essentially corre- 



sponds, in our language, to taking into account only the wave function corrections of 
eqs. ([25|j26|) . Indeed for m s of the order of 200 MeV the second order wave function 
correction almost exactly coincides with the exact result of the Yabu-Ando method. 
This is illustrated in Fig.4. However consistency requires to include the kinematical 



and dynamical contributions of eqs.(|22"D in the same order of m s . Whereas the kine- 
matical corrections are always small the dynamical ones are by no means negligible. 
This is explicitly shown in Tab. [Til], where all contributions to constants A ... H are 
displayed for M = 423 MeV and m s =180 MeV. Note that 0{m1) corrections to F are 
negligible. 

The message here is clear: The wave function contributions (w.f.) are not the whole 

17 



story and one has to include consistently all corrections, i.e. in addition kinematical 
and dynamical contributions, in a given order of m B . 

The purpose of this section was twofold: First we have demonstrated the impor- 
tance of the 0{m1) corrections coming from the effective action as compared to the 
wave function corrections (the Yabu-Ando approach). Second we have used the mass 
splittings to fix the parameter of the model, namely the constituent mass M. Having 
done this we can proceed to the evaluation of the axial coupling constants. 



IV. AXIAL CURRENTS IN CHIRAL LIMIT 

In order to calculate observables like the axial vector currents, one has to consider 
the path integral expectation value of these operators. This can be also done within 
the formalism of the quark correlation functions, see e.g. Refs. [[Ryn|]. Here will use 
however the effective action approach and show how the time- ordering within a quark 
loop together with the collective quantization brings up the corrections linear in the 
rotational velocity Q |T6|-[T8|. The approach presented here will be different from the 
one of Ref. fLTf , where these new corrections were calculated from the unregularized 
expressions. 

One can express the axial vector current A® as a path integral expectation value 

<Afa)> = N J VqVq J VU(q-f^ 5 X a q) e ~ I 

- J VqVq j PC/expj- J d 4 x (C E - s q W I a q)} | ^ (31) 



S 



Ss(x) 

with the convention 7 = —274 and the following definitions for I a : 

2i 0) :/ = l, 9k 3) :I = \ 3 , <?i 8) :/ = A 8 . (32) 

As in Sect. Ill the quantization is performed by introducing the time-dependent SU(3) 
rotation matrix A(t). Then we obtain 



< A a Jx) >= — - 
^ 8s(x) 



J V~qDq JvUexp^-J d A x (jC e - s q 7 ^X a Aq) } | s=0 (33) 



with q = A{t)^q and q = qA(t), and the rotated lagrangian £ E 
C E = ~q (-ip + M U(x) + A ] mA - i^A^A) q. 



(34) 



Integrating over the quark fields and restricting the T>U integration to the SU(3) 



rotations gives for the space components of the current 
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< Af(x) >-- 



8s(x) 



VA(t) S P(to) log D[s] 



s=0, 



(35) 



where D[s] = d t + H + A^A - ij A A^mA + i s 747^75 A^ \ a A. 

Expression ([35|) needs now a careful explanation. This is due to the fact that we 
are not going to perform the path-integral over the rotational matrices but, instead, 
we will use the operator formalism after the quantization of the generalized SU(3) 
coordinates has been performed. Within the path integral the time-ordering of the 
operators after the T>A integration is given in a natural way (see e.g. Ref. |36]]). In 
the operator formalism we can make use of the trace properties only if we respect 
the time- ordering of the operators (intrinsic and collective) within the trace. This is 
denoted in short by the modified trace Sp( to ) , which has the usual properties except 
for the time component. Only after the explicit time ordering of the operators the 
rotational frequencies Q can be again considered as time-independent. 

One great advantage of retaining the trace in this form is the straightforward ap- 
plicability of the regularization procedure. Let us split the effective action into real 
and imaginary part: 



< AUx) > 



VA(t) 
VA(t) 
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8s{x / 
8 ' 1 



Re Sp(to) log-D + i Im Sp (to ) log£> 

D 



8s{x) 2 



Sp(to) log D ] D + Sp^ log 



(36) 
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It is important to consider s(x) as explicitly time-dependent (see App. B). Preserving 
vector gauge invariance [3J| by using the proper time regularization |3(J we regularize 



the real part of the effective action: 

Sp ( to) log D*D -> Sp (to ) / ^<f>(u) exp (-uD*D) , (37) 

where <f>(u) is given in Sect. II. Note that for symmetric contributions [171, i.e. when 
the index of a generator J a is such that it commutes with the .D^-function, the time- 
ordering has no influence. 

In the following the separation into valence and sea part is done by introducing 
D' = D — n, where \i is a chemical potential with < \l < E va \ ||12|| , such that 
S%* = S$ + Sfg with the definitions S$ = S cS [D'} - S cS [D] and S S J? = S cS [D}. The 
subtraction of possible vacuum contributions is implicitly understood. 

A. The lowest order contribution ~ 17° 

The axial vector coupling constants g%, defined as the corresponding formfactor in 
the limit q 2 = 0, can be calculated from eq. (|36|). Noting that the functional integral 
T>A over the rotation matrices can be replaced by an ordinary integration / over 
the collective coordinates one can define an operator g\ such that 

g a A = J dU < B{U) I g% I B(U) >, (38) 



where | B(£a) > is the baryon wave-function of eq. (|23|) . Therefore g% is obtained 
by expanding < Af(x) > in eq. (|36D in terms of the rotational frequency and strange 
mass but without performing the DA-integration. The rotational velocity is then 
replaced by the generators of fll8f) yielding after integration over 3-dimensional space 
the collective operator g\ in terms of J a and Wigner functions D ab . 
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Therefore the lowest order result in Q (i.e. Q°) comes from the proper-time reg- 
ularized real part of the EEA (fffi). One obtains for a = 3 and 8 (see App. B for 
details): 

g a A (Q°,m° s ) = M 3 D a3 , (39) 

where A^I a A = D^Xf, for a = 3 and 8, and A^I a A = 1 for a = 0. At this level g\ = 0. 
The quantity M 3 = M 3 val + M 3 sea which comes from the real part of the action is then 
given by: 

M 3jVal = N c < v | 707375A3 I v > (40) 

and 

^3, sea = — ^ z2 <n\ 707375A3 I n > sign(E n )nY;{E n ) , (41) 
where the regularization function reads [pql : 

V7T Jo VM -E^ 



The values of M 3 are displayed in Tab. [V 



B. Anomalous 1/A^ C corrections from the imaginary part of the EEA 

Taking 1/N C corrections (i.e. terms linear in the rotational velocity Q) into account, 
one can make the clear separation into local quantities which emerge from the imaginary 
part of the effective Euclidean action, and non-local quantities emerging from the real 
part due to explicit time-ordering of collective operators. It will turn out that the 



former quantities are related to Witten's anomalous axial current |53 |, whereas the 
latter ones have no counterpart in a mesonic effective theory. In a certain sense they 
renormalize the leading contribution of the axial current given by eq.(|3"9"D (compare 



with Refs. H,|4l|). 
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In the chiral limit the anomalous corrections linear in f2 can be written as: 



g% = -M bc i{D ab , Q C E } = -2M bc D ab Q c M (43) 

with 

M bc = -J- Yl < n I a ^ I m >< m | A c | n > K M (E n , E m ) (44) 

^ n,m 

and the cutoff independent regularization function 71m 

zip 7~n \ 1 sign(ff ra - fj) - sign(ff m - /i) 
TZ M {E n ,E m ) = _ . (45) 

As noted above the chemical potential fi is chosen in such a way, that it always lies be- 
tween the valence level and the positive continuum of states. In this way the quantities 
M bc correspond to the full baryon number one contribution and therefore contain the 
sum of the valence and the sea part. Additionally we define for later use M 8a = ^/3M 8a 
and M a8 = \/3M a8 . The only non- vanishing contributions in eq. (|43| ) are: M 83 and 
M 44 = M 55 = — M 66 = — M 77 . Using the symmetries of the hedgehog states one can 
write for a = 3 and 8: 

g% = -—^D a8 n 3 - AM u d 3bb D ab Q b , (46) 

where the sum over b = 4 ... 7 is understood. The values of the coefficients entering 
eq.p6|) are displayed in Tabs. [TV] and [VI]. It is clear from the form of eq.fl4"6|) that the 
anomalous corrections linear in Q vanish in the SU(2) case. 

C. Normal 1/N C corrections from the real part of the EEA 

Apart from the anomalous contributions of the preceding subsection, which have 
their counterparts within the SU(3) Skyrme model f42| , |2"4|1 , there exist corrections linear 
in Q, which come from the real, non-anomalous part of the EEA and which are due 
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to some explicit time-ordering of operators. They were recently discussed within the 
SU(2) |l6l , |T8l1 and SU(3) [17] version of present model. As advertised in the beginning 
of this section these terms emerge because the operators J a and D ao in the space of the 
collective coordinates are in general time-dependent and have to be time-ordered, if 
one makes use of the operator formalism . This is described at length in Appendix 
B and as a result the contribution to g A from these terms can be summarized as: 

9l = ^ [Wk, D ab ] E< Jl U>><m a-X \ n > TZ Q (E n , E m ) (47) 

^ m,n 

where the rather complicated regularization function TZq is given by: 

K (E E ) = f 1 da aE n-{ l -<x) E ™ exp (-[ag + (1 - a)Ej}/ 'Kj) 
Q{ m> Jo 2vr J a (i- a ) 1 aEl + {l-a)El * 1 ] 



Here the proper-time w-integration for our step-like functions (j)(u) has been already 
performed (see App.B for a general expression). In the limit Aj — ► oo eq. (^) im- 
mediately reduces to eq. ( B13|) of App. B and coincides therefore with our former 



prescription in Ref. |l~Tfl. However, as we will see later, with the regularization prop- 



erly taken into account, the physical values will come out much better. 

Using the quantization condition for Q and making use of the commutator 



if c bdD a d |43fl eq. fl47|) can be written as 



~ a ifcdbDad 

9a = t " bc 



1 cc 



2iQn + . ^ 



h h 

where the quantities Qb c coming from the real part of the EEA are given by Qb c = 
Qfec.vai + Qfec.sea- Explicitly the valence part reads: 

Qfecvai = — 2^ B B signK (50) 

and the sea part: 
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Qbcna = ~Y £ < n I 7o7375A b | m >< m | A c | n > TZ Q (E n , E m ). (51) 



The numerical values for the Qb c can be found in Tab. |V|. 

One should note however that the appearance of Qb c is due to the fact that the 
operators in the EEA are explicitly time-ordered and that in addition the matrix 
elements < n | 707375 | m >< m | A c | n > are antisymmetric with respect to the 
interchange of m and n. 

The valence contribution Qb c ,vai differs from the formula given in Ref. [pj, where the 
existence of such corrections was claimed for the first time. The correct path-integral 
formula is given by eq. ([50]) and eq. (|5T|). Numerically however the difference between 
our expression for Qb c ,vai and the expression of Ref. |16| is quite small. Note also that 
in Ref. |16| the sea contribution to Qbc was erroneously claimed to be identically zero. 
Again numerically Q bC)S e& is rather small. 

Putting all these corrections together one obtains: 

2iQi2 2iQ 45 ~ 



g a A (Q + Q\m° s ) 



h l. 



2 



D 



«,3 



2M 83 f M 44 f 

D a8 J 3 - 4—— d m D ab J b (52) 



(b runs over 4 ... 7). Note that all the quantities M3, Q bc , Mb c and also Ji^ are of the 
order 0(N C ), such that the Qb c -terms in the brackets indeed correspond to 1/N C cor- 
rections to the lowest order result. In other words as far as one neglects the anomalous, 
purely SU(3) contribution in eq. ([52]), the ratio of different g^s has no 1/N C correction. 

D. The anomalous singlet axial current 



The singlet axial vector current was already given in [53] and it gets only anomalous 
contribution linear in Q: 

-0 /nl 0\ 2M 83 f 

g A (Q ,m s ) = — J 3 . (53) 

i\ 
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Note that eq. ( |53[) given here in the context of SU(3) coincides exactly with the SU(2) 
result of Ref. This is because only spin eigenvalues (J3) enter here, whereas the 
other g^s always contain D-functions, whose matrix elements depend crucially on the 
SU(iVfl avor ) algebra used. 



E. The axial currents in the leading order LWLA 

For large size of the soliton one can approximate the different contributions in eq. 
|)2|) by the gradient expansion (or long wave-length approximation (LWLA)) ||45|| . The 



lowest order result for SU(2) is given in Ref. |22] and it coincides with the expressions 
from the Skyrme model. Terms linear in Q can be also gradient-expanded. In SU(3) 
one gets the anomalous contribution coinciding with the Wess-Zumino-Witten term in 
the Skyrme model. Using the results of App. C one can also calculate the LWLA of 
the time- ordered termsQ. Altogether one obtains: 

sin 26 N 



9l = Jdrr 2 le' + 



8?r ,» M M ' 

— f H 1 

3 J * Ah 8/ 2 . 



D 



«3 



4 N r 

drrO 1 sin 3 9{r) d 3bb D ab J b . (54) 

lo 07T J 



4 N. 
h 

It is clear from eq. ([54]) that the time-ordered contributions (i.e. Q bc ) lead to the 
renormalization of in the sense that they are also proportional to D a3 . In SU(2) 
where the second line of eq. ([54]) vanishes the ratio of ^'s for different baryons has 
no 1/N C correction. This was also found by Dashen and Manohar from large N c QCD 
pO| , ^T| . Furthermore the SU(2) result resembles very much the old non-relativistic 
quark model prediction for the 1/N C correction which is given by: g& = N c /3 + 2/3. 
Using the value h = 1.156 fm for M = 423 MeV from Tab. |, the SU(2) part of eq. 



3 The moments Q bc are evaluated here in the infinite cutoff limit in order to get a simple 
and cutoff independent result. 
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(]54T) gets approximately 50% correction from the 1/N C term. In SU(3) however, there 
are additional corrections from the third term in the brackets of eq. ( |54"D and from the 
anomalous terms, such that the simple rescaling factor does not exist any more. 



V. MASS CORRECTIONS FOR G A 

In this Section we will evaluate the symmetry breaking corrections to the axial 
currents due to the non-vanishing strange quark mass. These arise from the term 
A^mA = /io — [i%\ a Dsa- In the linear order in m s and in the zeroth order in Q neither 
contributions from the imaginary part nor from the explicit time-ordering (because D- 
functions always commute with each other) exists. Therefore entire symmetry breaking 
contribution comes from the real part of the EEA. Performing the expansion of the 
real part of the EEA in m s one gets: 



4m s „ „ „ , 4t7\ it 

"83 



gli^M) = --^R 38 D a3 (l - D 88 ) + ^R 83 D a8 D 
8m. 



d 3 bbRuD a bDsb (55) 



V3 

with b = 4 ... 7 The proper time regularized quantities R bc = Rbc,va\ + Rbc,sea are given 



by: 



and 



N c ^ < n | a 3 \ b | v >< v | A c 7o | n > 

RboM = ~7T 2^ f _ f ( 56 ) 



Rbc,se* = -j- < 71 I °"3A 6 | m >< m I A c 7o | n > Kp(E n , E m ). (57) 

n,m 



with 



n,(E n ,E m ) = -^= ^m\ n F f ] (58) 

2^1T Jo y/t E n - E m 

For future use we also define R$ 3 = V^Rss and R 3 % = \/3R 3 8- Note that 1lp(E n , E m ) 
is different from the regularization functions TZj(E n , E m ) and TZq(E u , E m ). The origin 
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of this difference, which however survives only in the finite cutoff case, is the different 
hermiticity behavior of the current and the mass term on the one hand and the Coriolis 



term zQe in eq. (|35|) on the other hand. The latter one turns out to be antihermitian 
in Euclidean space, whereas the former ones are hermitian. Because the proper time 
regularization rests on building -D e -De from the very beginning, different signs emerge 
and lead to different regularization functions. Their substantial different behavior can 



be seen in Fig. f of Rcf. [TT] and Fig. 5 of the present paper. We list the new coefficients 
R a b in Tab. [VT] and Tab. Vll . 



Apart from these terms from the action we have in addition the m s terms from 
the quantization condition ([181) . Including these mass corrections and the ones from 
the expansion of the effective action we obtain up to the linear order in the symmetry 
breaking and the rotational frequency (for a = 3 and 8): 



91 



2iQ 12 2iQ i5 



4M 44 j 

L> a3 r a 3bb- u ab<Jb 



2M 83 4m. Id j 
4m v 



V3 h 



V3 
4m s 

'73 



R^Da&Ds'i + 



5f7l 



—7=-(Rm — M44——)d S bbDabD8b 

V3 h 



R 38 D a3 (l-D 88 ), 



(59) 



where, as usually, the index b in g^;, runs over 4 ... 7. 

The quantities R$3 and M§3 are already known from the expression of the flavor 
singlet axial constant [ PJ . We found thereQ in the same order of the expansion: 



9l 



2M 



S3 



J 3 



4m s 



-D 



83 



The main difference between eq. (1591) and eq 



K x — - 

— — M 83 — i? 83 



(60) 



is that the lowest order term for 



Comparison with Ref. |44] can be done by identifying = (3\ and R$z = @2- Note also 



that the sign in Eq.(5) in Ref. [ 44 1 is misprinted 
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g A is purely anomalous, whereas the corresponding term for the g\ is non-anomalous. 
We come back to this point, when we make the comparison with the Skyrme model in 
Sect. VI. 



VI. NUMERICAL RESULTS FOR AXIAL CURRENTS 



The three different measurements of the spin asymmetry in the polarized lepton- 



nucleon deep inelastic scattering |46]- |I9| have been recently reexamined by Ellis and 



Karliner ||50|| . The message of their work is that whereas the Bjorken sum rule [pTJ is 



in agreement with the data, the Ellis- Jaffe sum rule |52|J53[| is violated and the results 
read finally: 

gf' exp = 0.24 ±0.09, gf' exp = 0.35 ± 0.04 and gf' exp = 1.25. (61) 

In this section we discuss our numerical results for the three axial decay constants 
9a\ 9a^ an d 9a^ including the strange mass corrections. They are summarized in 
Tab. [VIII| , Tab. [IX| and Tab. [X]. Our final values for a constituent quark mass 



M = 420MeV are given by 

= 0.37, #| 8) = 0.31 and gf = 1.38. (62) 

(3) 



In Tab. | V 1 1 1| the difference between SU(2) and SU(3) results for g A ' can be 



seen in each order of the 1/N C expansion. Obviously the lowest order contribution 
(Q°) in SU(3) is significantly smaller than in SU(2) due to the fact that the SU(3) 
expectation value of the corresponding D-function D 33 is only 70% of the SU(2) value. 
The anomalous contribution of Eq.([5J) linear in Q which is non-zero only in the SU(3) 
case acts as a substitute for this group-theoretical reduction. Indeed, it leads to an 
almost exact readjustment of the SU(3) value to the SU(2) one. For our preferred value 
of M = 423 MeV from the hyperon spectra and m s = 180 MeV it pushes the leading 
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order SU(3) result up to g A ~ 0.84 . These two values of the model parameters M 
and m s will be used in the following discussion of the numerical results. Due to the 
presence of the quantities Qb c from the explicit time-ordering the SU(2), as well as the 
SU(3) results, have corrections linear in the rotational velocity. These conceptually 
new terms have no counterparts in the ordinary Skyrme model. Similarly to the old 
non-relativistic quark model estimates of the 1/N C correction, i.e. g A = N c /3 + 2/3, 
these new terms turn out to be of the order of 50% of the leading term. For SU(2) 
they push g A from 0.84 to 1.15 and in SU(3) they give the final value of 1.31. Note 
that the latter value is obtained with regularized time-ordered quantities Q^ c . In the 
present work we have derived Qb c from the proper-time regularized real part of the 
EEA. In Ref. ||17|| , where the regularization was neglected as a first approximation f], 



the sea part of the quantities Qb c made a ~ 30% contribution to the total value of 
the Qbc- Although the Qb c are finite it should be regularized since it is a quantity 
coming from the real part of the EEA. This is done in the present paper and as a 
result with the inclusion of a regularization function the sea contribution is less than 
3% (for M = 423MeV, see Tab. 0). 

Various contributions from the strange quark mass (kinematical, dynamical and 
wave function) increase the value of g A of about ~ 5% up to g A = 1.38, such that 
the experimental value g <y f >,exp = 1.25 is overestimated only by ~ 10%. 

It has to be stressed that this is in contrast to all calculations within the purely 
pseudoscalar [M] or pseudoscalar and vector Skyrme model [2"4"|,l"2"||, in which g A is 



underestimated by ~ 30%. That this significant difference is due to the presence of 
the new terms from the time-ordering of the functional trace is most clearly evident 



5 Note that for this purpouse the functional trace Sp had to be changed to Sp( to ) which 
has different properties in time indices. 
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from Tab. [X|. There the flavour contributions to the axial current are given for the 
Skyrme model and the NJL model with and without the time-ordered (T) corrections. 
Without the new corrections the NJL model resembles very much the numerical results 
of the SU(3) Skyrme model with vector mesons. This was already noted at the level 
of the collective hamiltonian for the mass splittings in Ref. [[HJ and here can be seen 
numerically for the axial currents with a high accuracy. 

Apart from g A we list also the values for g A \ partially given already in Ref. [P], 
and g A . The spin of the proton, that is carried by the quarks, and is equal to the 
matrix element of the flavour singlet axial vector current, has no contribution in the 
order Q°, but gets a first non-vanishing contribution in the linear order of Q. This, 



as can be seen from Tab. [IX|, is also the dominating contribution which gets only 
a very small strange mass correction from the kinematical or dynamical parts. The 
corrections from the wave-function vanish in the linear order of m s , because the spin 
operator is diagonal in space of the higher SU(3) representations. For M = 423 MeV, 
the theoretical value of g A ~ 0.37 is a little bit above the experimental error bars, 
nevertheless one has to keep in mind that the analysis of the experimental data is still 
under debate |50] . 

(8) 

Experimental extraction of g A from the hyperon semi-leptonic decays depends on 
how the strange quark mass corrections to the SU(3) symmetric result are taken into 

account. Therefore the experimental error bars on this quantity may be at present 

(s) 

too small. In the present calculation we obtain g A = 0.31 to be compared with the 



'experimental' number of |5(| g < fi >,exp = 0.35 ± 0.04. 

So from our calculations one can conclude that for the 'fixed' mass of M = 423 MeV 
and m s = 180 MeV all three axial vector coupling constants are quite close to the 



experimental values of [|50|. From Fig. 6 it can be seen that for larger mass of M ~ 
550 MeV, g A and g A almost coincide with the experimental values, whereas g A \ 
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having relative large negative strange quark mass correction, deviates from the central 
value 0.35. However, ignoring the large m s corrections (as it is usually assumed in the 
analysis of the hyperon semileptonic decays @]), even the value of coincides with 
the experimental value of Ref. |50| . 



In our view there is no reason to talk of a spin crisis of the nucleon. In the present 
approach the proton is treated as a many-body system of valence and sea quarks and 
the data of the EMC-experiments are basically reproduced. 

VII. COMPARISON WITH THE SKYRME MODEL 

Now we want to compare our results with the Skyrme model, which can be regarded 



as a large constituent quark mass limit of the NJL model |L2|, or - equivalently 



the large soliton size limit. It was already mentioned in Ref. p5] that the formula 
for the mass splitting of the present SU(3) NJL model is much more subtle than 
the corresponding expression for the pseudoscalar Skyrme model. This is because 
certain important anomalous terms dominated by the valence quarks are missing in 
the effective action of the Skyrme model. This has been cured in the Skyrme model 



only by introducing vector mesons Therefore we will focus here on the Skyrme 

model, in which vector mesons and in addition kaon fluctuations and the gauged Wess- 
Zumino term are added. Then the collective operator has the following structure 



9 A — a \D 33 + a>2d3aaD3 a R a + a 3 D 38 + a/id 3aa D 3a D 8a 

+a 5 D 33 {l - D S8 ) + a 6 D 38 D S3 (63) 

which corresponds effectively to the expression for the NJL model. Although the 
origin of the various coefficients are quite different in the NJL and Skyrme model, 
both approaches give effectively the same operator structure for g^. 

However one should stress here two important differences: first of all the new 
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corrections linear in which arise due to the time- ordering within the fermion loop 
vanish in the Skyrme model identically. The Skyrme model is based upon a local 
lagrangian density which, apart from the Wess-Zumino term, is even in time derivatives 
and therefore the spatial components of the axial currents are also even not allowing 
for terms linear terms in Q. Second, even if one restricts oneself to the terms not 
including the corrections due to the time-odering (local limit) the contribution of the 
valence quarks in the present model makes our results qualitatively very different from 
the ones of the Skyrme model. The coefficient a 2 e.g. is in the Skyrme model with 
purely pseudoscalar mesons dominated by the induced kaon fluctuations |23|, which are 
neglected in the present NJL model. If the vector mesons are included in the Skyrme 
approach the situation does not change qualitatively [p4| . In the Skyrme model a 2 



gives only a 10% contribution to g& [ 23j| , whereas the dominating valence contribution 



to M44 (compare Tab. |V|) in the NJL-model gives almost a 30% contribution to gx if 
the terms due to the time- ordering are neglected. The fact that the total values for g^ 
in the local limit of the present approach and in the vector meson pseudoscalar Skyrme 



model in Ref. [23] roughly coincide hinges also on the rescaling procedure for the 
parameter e in Ref. p3"| , which tends to increase g^. The Wess-Zumino contributions 
to (?a in the Skyrme approach play a minor role, in the NJL model the anomalous part 
of the action, containing the WZ term in lowest order of the gradient expansion, gives 
~ 1/4 of the total amount. 

VIII. SUMMARY AND DISCUSSION 

In this paper we have investigated mass corrections to the hyperon masses and to 
the axial currents of the semibosonized SU(3) Nambu-Jona-Lasinio model. Moreover 
we have extended our recent analysis of the corrections to the axial currents which 



appear due to the time-ordering of the quark loop and semiclassical quantization [[17 1 
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to the case of the regularized effective action. In the semibosonized NJL model baryons 
are understood as solitonic solutions of the classical equations of motion. However 
the solitons do not carry proper quantum numbers and the semiclassical quantization 
procedure has to be applied in order to describe the mass splittings within the strange 
baryon multiplets. This treatment is based on introducing time-dependent rotations 
in the direction of the zero modes followed by a canonical quantization of the collective 
coordinates of the rotation matrix. Since these zero modes contribute significantly to 



the mass splittings [20], it was a challenging task to look at the axial currents which 
can be related to the recent measurements of the spin structure functions ^D,E7|,ff^,HH] . 

First the model parameters were fixed by looking at the hyperon mass splittings up 
to the terms quadratic in m s . These are reproduced with unexpectedly high accuracy 
and point towards a constituent quark mass of M ~ 420 MeV. We have also explicitly 
shown that the wave function corrections and the corrections due to the expansion of 
the effective action are comparable and therefore it is somehow inconsistent to perform 
only the Yabu-Ando diagonalization of the first order hamiltonian. 

Second we considered the axial vector currents with the inclusion of the linear 
corrections in the rotational velocity. The new contributions which appear due to 
the time-ordering of the quark loop and semiclassical quantization have been shown 
explicitly to come from the real part of the effective action. If the proper-time regular- 
ization is implemented they are dominated by the contribution of the valence quarks. 
In the SU(3) model there are also other contributions linear in Q which come from the 
imaginary part of the effective action and as such do not require regularization. They 
are also dominated by the valence contribution. This concerns the leading term of g ^ , 
which vanishes in the pure pseudoscalar Skyrme model, whereas it is non-vanishing 
(however small) in the present model in rough agreement with experiment. 

The expression for g A has a ~ 25% rotational contribution from the imaginary 
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part of the effective Euclidean action, which vanishes in the SU(2) case and which 
can be related to Witten's formula for the axial vector current from the Wess-Zumino 
effective action. Moreover it has a ~ 30% contribution due to the explicit time- ordering 
{Qbc) of the collective operators. These terms are not present in local theories like the 
Skyrme model. In the present model they are entirely due to the non-locality of the 
fermion determinant. Performing the gradient expansion of these quantities, it can be 
shown that these terms have the same mesonic structure as the lowest order term, such 
that one can define a profile independent renormalization of the lowest order g A . This 
is similar to recent findings of Dashen and Manohar JIT] within large- iV c QCD and to 
the old non-relativistic quark model result of = (N c + 2)/3. 

Quantitatively despite the fact that the lowest order SU(3) result is reduced by a 
group theoretical factor of 0.7 with respect to the SU(2) case the new time-ordering 
and anomalous contributions push the total value of g A upwards. 

Next we have considered corrections linear in the strange quark mass. They are 
consistently derived from the expansion of the effective action, the quantization condi- 
tion as well as from the higher representations of the wave function in the spirit of the 
Yabu-Ando diagonalization. However the effect on g A is not large and finally one ends 
up with g A = 1.38 for M = 423 MeV, which is only ~ 10% above the experimental 
value of g^' exp = 1.25. Here it should be stressed, that such nice agreement was never 
obtained within the pseudoscalar or pseudoscalar vector Skyrme model |2"3"| , |2"4]] . This 
qualitative and quantitative difference comes from the new 1/N C corrections of the 
time-ordering from the non-local structure of the present NJL model. 

The g^ exists already in SU(2) and the only effect in SU(3) is a small shift due 
to the finite symmetry breaking m s . This is in contrast to g A , which vanishes in 
the SU(2) case, and which in SU(3) gets entire contribution from the rotation and 
from the strange quark mass. In the chiral limit the values for g A and g A are quite 



34 



close to each other, as suggested in [55|], however the strange mass corrections reduce 
the value of by ~ 25%, whereas the explicit symmetry breaking has almost no 
influence on This holds at least if we take all linear m s corrections into account 
and even the m^-corrections, which can be calculated in this framework from the 
non-symmetric wave- functions [21]. The final values g^' ~ 0.37 and g^ ~ 0.31 for 
M = 423 MeV and m s = 180 MeV are to be compared with the experimental data 
extracted from the recent Ellis and Karliner analysis |5(J, i.e. g^' exp — 0.24 ± 0.09 
and g^' exp — 0.35 ± 0.04. Apparently the theoretical values are only slightly outside. 

Qualitatively the following can be said: g^, which represents the part of the proton 
spin carried by the quarks, gets a non-vanishing expectation value entirely due to the 
anomalous part of the EEA. In a constitituent quark model, when the total spin of the 
proton is entirely carried by three quarks, this g^ equals one. The present model gives 
values close to experiment since the proton is treated entirely as many body system 
rotating in spin-isospin space. Hence one cannot attribute the spin of the proton 
to single elementary particles but only to the whole system. Hence the fact that in 
agreement with the experiments only a fraction of about 20 — 30% of the nucleon spin 
is carried by the quark-spins is not surprising at all in the present model. One should 
note that in the Skyrme models a non-vanishing value of g^ can be obtained only by 



adding vector mesons for the anomalous part E 

Altogether the picture which emerges is quite satisfactory. Mass splittings are 
accurately reproduced and axial currents are in good agreement with the experimental 
data if rotational 1/N C corrections are taken into account. In particular the spin of the 
proton originates in this model to about 35% from the spin of the quarks, a number 
being in reasonable agreement to the world data reported by SMC j|8|]. Together with 



the numerical results for the Gottfried sum |57|] the model provides a good reproduction 
of experimenta data so far. 
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On purely theoretical side the presence of the new terms linear in f2 calculated in 
this paper poses a serious problem to effective meson theories like the Skyrme model, 
where such terms vanish identically. Another theoretical question which deserves a 
comment is the convergence of the expansion in Q. The large size of the corrections 
calculated in this paper, although expected from the quark model calculations, poses 
a problem of reliability of the numerical results. One way to tackle this problem would 
be to calculate the Q 2 corrections to the axial currents. This is the highest power of 
Q one can think of, since the collective hamiltonian itself is of that order. Despite 
the technical difficulties in calculating these terms the preliminary estimates indicate 
that they are not negligible^]. Therefore one has to conclude that the expansion in 
Q is slowly convergent. Moreover the formalism of the collective quantization has to 
be revised if one wants to include terms higher than Q 2 . In addition despite the fact 
that mass splittings are well reproduced the absolute energies provide still some prob- 



lems which are associated with zero-point corrections p8| , p0| and boson fluctuations. 
These questions are certainly beyond the scope of this paper, where we had to content 
ourselves with the linear corrections alone. 
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TABLES 



TABLE I. Moments of inertia for different constituent masses 



M 


£ [SU(2)] 


h 


h 


K x 


K 2 


N 




N 2 


[MeV] 


[MeV] 


[fm] 


[fm] 


[fm] 


[fm] 


[fm] 


[fm] 


[fm] 


363. 


60.32 


1.512 


0.720 


0.606 


0.372 


0.765 


0.647 


0.496 


395. 


58.14 


1.285 


0.618 


0.438 


0.290 


0.704 


0.500 


0.408 


419. 


56.14 


1.178 


0.569 


0.369 


0.255 


0.668 


0.438 


0.370 


423. 


55.52 


1.156 


0.560 


0.357 


0.250 


0.658 


0.426 


0.362 


465. 


51.86 


1.015 


0.496 


0.276 


0.210 


0.599 


0.349 


0.311 



TABLE II. Coefficients of Eqs.(22) as functions of constituent mass M for m s = 180 MeV 



M 


A 


F 


D 


G 


B 


C 


H 


363. 


489.04 


427.02 


94.67 


1.96 


450.89 


510.24 


6.70 


395. 


481.32 


399.60 


97.41 


1.70 


449.48 


457.80 


5.56 


419. 


469.60 


381.66 


97.17 


1.51 


441.70 


427.02 


4.86 


423. 


465.57 


377.47 


96.79 


1.48 


438.60 


420.32 


4.70 


465. 


441.39 


349.75 


94.91 


1.25 


419.82 


376.11 


3.72 



38 



TABLE III. Different contributions to coefficients of Eqs.(22) for a constituent mass 
M = 423 MeV and for m s = 180 MeV 



0(m s ) 








0(ml) 




total 


exp. 




kin. 




dyn. 


w.f. 


total 






A 546.10 


10.94 




64.64 


-0.15 


-53.85 


492.25 


— 


B 546.10 


10.85 




64.55 


-53.81 


-107.50 


438.60 




F 381.20 


1.18 




27.67 


22.76 


-3.73 


377.47 


379.00 


D 120.76 


-0.02 




11.78 


-0.07 


-11.87 


108.89 


86.00 


C 348.16 


1.20 




19.97 


90.92 


72.15 


420.32 


422.00 


G 0.00 


0.61 




-0.66 


1.53 


1.48 


1.48 


6.75 


H 0.00 


-0.29 




0.41 


4.57 


4.70 


4.70 


2.20 


TABLE IV. 


Quantities M3 and M. 


44 for the SU(3) Nambu-Jona-Lasinio model 


in de- 


pendence on the constituent quark mass M 












M 3 








M 44 [fm] 




M [MeV] 


val 


sea 




tot 


val 


sea 


tot 


363 


-2.293 


-0.468 




-2.761 


-0.414 


-0.011 


-0.425 


395 


-2.235 


-0.384 




-2.619 


-0.326 


-0.012 


-0.338 


419 


-2.209 


-0.316 




-2.525 


-0.288 


-0.012 


-0.301 


423 


-2.205 


-0.307 




-2.511 


-0.283 


-0.013 


-0.295 


465 


-2.173 


-0.203 




-2.377 


-0.238 


-0.014 


-0.251 
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TABLE V. 


Quantiyies Q12 = 


-2iQ_ + 


and Q45 = 


for the 


SU(3) 


Nambu-Jona-Lasinio model in dependence on the constituent quark 


mass M 








Q-+ M 






Q45 [fm] 




M [MeV] 


val 


sea 


tot 


val 


sea 


tot 


363 


.408 


.025 


.433 


-.410 


-.023 


-.433 


395 


.317 


.021 


.339 


-.318 


-.020 


-.338 


419 


.279 


.019 


.298 


-.279 


-.018 


-.297 


423 


.272 


.019 


.291 


-.273 


-.017 


-.z90 


465 


.226 


.016 


.242 


-.226 


-.015 


-.241 


TABLE VI. 


Quantities M§3 and i?§3 for the SU(3) Nambu-Jona-Lasinio model 


in de- 


pendence on the constituent quark mass 


M 












M 83 [fm] 






R 83 [fm] 




M [MeV] 


val 


sea 


tot 


val 


sea 


tot 


363 


-0.683 


-0.015 


-0.699 


-0.293 


-0.094 


-u.oo / 


395 


-0.500 


-0.016 


-0.516 


-0.150 


-0.090 


-0.240 


419 


-0.422 


-0.016 


-0.438 


-0.095 


-0.091 


-0.186 


423 


-0.409 


-0.016 


-0.425 


-0.086 


-0.091 


-0.177 


465 


-0.316 


-0.017 


-0.333 


-0.025 


-0.090 


-0.115 
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TABLE VII. Quantities R44 and R38 for the SU(3) Nambu-Jona-Lasinio model in de- 
pendence on the constituent quark mass M 



M [MeV] 




i? 44 [fm] 






R 38 [fm] 




val 


sea 


tot 


val 


sea 


tot 


363 


-0.253 


-0.023 


-0.277 


0.074 


-0.049 


0.024 


395 


-0.178 


-0.029 


-0.207 


0.083 


-0.066 


0.017 


419 


-0.148 


-0.030 


-0.179 


0.086 


-0.073 


0.012 


423 


-0.144 


-0.031 


-0.175 


0.087 


-0.075 


0.012 


465 


-0.111 


-0.033 


-0.143 


0.090 


-0.083 


0.006 



TABLE VIII. The axial vector coupling constant g A for the SU(3) Nambu-Jona-Lasinio 
model in dependence on the constituent quark mass M. The strange current quark mass is 
chosen as m s = 180 MeV. The final model predictions are given by g^^ 1 ) in SU(2) and 
g^\Q 1 ,m s ) in SU(3). The experimental value is given by g^' exp = 1.25. 



M [MeV] 


SU(2) 






SU(3) 




gfis?) 9fm 


gf{UL\ml) 


«7i 3) (n\m°) 




363 


0.920 


1.302 


0.644 


1.482 


1.603 


395 


0.873 


1.224 


0.611 


1.381 


1.473 


419 


0.841 


1.179 


0.589 


1.328 


1.407 


423 


0.837 


1.173 


0.585 


1.314 


1.380 


465 


0.792 


1.109 


0.554 


1.250 


1.308 
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TABLE IX. The axial vector coupling constant g A and g A for the SU(3) 
Nambu-Jona-Lasinio model in dependence of the constituent quark mass M. The strange 
current quark mass is chosen as m s = 180 MeV. The final model predictions are given by 
p^(n 1 ,m s ) and g^(n 1 ,m s ). The experimental values are given by g^' exp = 0.35 ± 0.04 



and g^>' exp = 0.24 ± 0.09 (Ellis and Karliner [50]). 

M [MeV] gj 8) (O ,m s °) ffg(gVrg fff^X) gfgVj) 

363 0.159 0.443 0.328 0.462 0.475 

395 0.151 0.408 0.316 0.401 0.409 

419 0.145 0.389 0.309 0.371 0.377 

423 0.144 0.385 0.308 0.364 0.370 

465 0.137 0.363 0.301 0.328 0.331 



TABLE X. Various contributions to the axial vector current of the proton in terms of 
u,d and s quarks. A comparison is made between the Skyrme model with vector mesons [24] 
(Skyrme,vector), the NJL model without (NJL,scalar) and with the time- ordered corrections 
of this paper (NJL, scalar, T). In the last column experimental values from Ellis and Karliner 
[50] are given. 

Skyrme (vector) NJL (scalar) NJL (scalar, T) 'experiment' 

An 0.63 0.64 0.902 0.806 

Ad -0.31 -0.24 -0.478 -0.444 

As -0.03 -0.02 -0.054 -0.120 
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TABLE XI. Matrix elemnets of the operators for g\ in the proton state with spin up, 
where the index i is always running from 1 to 3 and b from 4 to 7. 



^33 


^38 


^88 


dmDzbJb 


d3bbD$bJb 


d3bbD3bD$b 


d3bbDsbD 8 b 


-7/30 


V3/30 


3/10 


7/60 


1/(20^3) 


-11/(90^3) 


1/30 


£>38£>83 


-C33-C88 


^83 


-C83-C88 




D 3t D 8i 




-1/45 


-4/45 


-^3/30 





7/20 


^3/45 


-^3/45 
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APPENDIX A: SEMICLASSICAL QUANTIZATION 

In this Appendix [] we collect all essential ingridients of the semiclassical quanti- 
zation scheme. For the purpose of simplicity we will consider only L without mass 
corrections. The generalization to include m s terms is straightforward. To this end 
one defines curvlinear angular velocities: 

A* A = '-ttiXt = l -9 k e kl {9)K (Al) 

where the vielbeins fulfil: e k i e k j = 8ij and ejf. = 8ij (the bar denotes the inverse 
vielbein). The generalized momenta are defined as: 

8Lq 

where c = N C J \[Yl. Here one postulates canonical commutation rules: 

[rTi, 9j] = -iSij (A3) 

and then the differencial representation for the operator fr is, as usual: 7Tj = —id/d6{ = 
-idi. 

It is convenient to define the quantities J,: 

Jj — eijiTi = VL m I mj - c S 8 j. (A4) 

The normal quantization procedure is here slightly subtle, since the tensor I is not 
invertible and hence Jg = — c is a constraint rather than the dynamical variable. To 
circumvent this difficulty let us introduce a small moment of inertia Jgg = e and then 
take e — > 0. Then one can proceed in a normal way, inverting Eq . (|A2|) for velocities 
and calculating the hamiltonian: 



7 We are greatful to P. Pobylitsa for discussion and for making clear to us many subtleties 
of the group theory 
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# (0) = M cl + -L £ J} + ±- 2 J? + ^( J 8 + c) 2 . (A5) 

i=1 i=4 ^ 

In the limit £ — > one recovers the constraint. 

Upon quantization Jj are promoted to the operators: Jj = —i e^di. It is straigth- 
forward, but tedious to calculate the commutation relations for Jj's. To this end one 
has to use the following identities: 

CjmdjGin Gj n 0j6i m &ikfkmn- (A6) 

Finally one arrives at: 

[J» Jj} = ifijkJk- (A7) 

The next question is how to interpret operators Jj. There are two global symmetries 
of the rotation matrix A which do not change the Lagrangian L : 

left : A -> e'ffi A and right : A -> A e'^ 1 . (A8) 

In the case of the right multiplication only Ai )2 ,3 enter, since due to the form of the 
hedgehog Ansatz the right multiplication has to commute with the right multiplication 
by A 8 . Then it becomes clear that the right symmetry corresponds to rotations, since 
it can be undone by the rotation of vector n entering the hedgehog Ansatz. On the 
contrary the left multiplication is just a global symmetry of the Lagrangian and it can 
be interpreted as a flavor symmetry. To calculate the generators of these symmetries 
one uses the Noether construction promoting ip to a time dependent quantity. Then, 
for infinitesimal transformations, one gets that: 

left : Qi — > Qi + (fijDji and right : Sl^ — > Sl^ — 

and subsequently: 
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d d 
left : - 7— L = -DijJj = Ti and right : - -7— 

ocpi d(fi 



L — Jj. 



(A9) 



Here D, tj = D$ = 1/2 Tr(A^XiA\j) are the Wigner matrices in the adjoint represen- 
tation of the SU(3) group. We have already shown that upon the quantization Jj are 
promoted to the SU(3) generators, however only Ji,2,3 correspond to the symmetry 
generators, namely to spin. To evaluate the commutation relations for Tj one has to 
convince onself that: 



Then it follows that: 



[Ti, Tj] — ifijkTk- 



(A10) 

Making use of the first of Eqs. ( |A9| ) and of the ortogonality relation for the D functions: 
D ac D bc = D ca D cb = S a b one can rewrite as in eq. fl2~T|). 
Baryon wave functions should fulfil the following identities: 

exp(iuj a f a ) tfj] (A) = tj; (V^ A "/ 2 v4) , [exp(ioo a J a ) ip] (A) = tjj (Ae iu}aXa/2 ) . (All) 



The phase convention in Eq. (|All|) is chosen in such a way that the ^'s are faithful, 



i.e.: 



exp( 



(too a 2) f a ) exp(^i 1) f a ) (A) = $ ( { e lJ ' )Xa/2 e iuJ «' Xa/2 \ A 



{iuj { a ^J a ) exp(z4 1) J a ) 4j {A) = tfj(Ae iJ ' )Xa/2 e iuJ »' Xa/2 



,(!)• 



exp( 

The problem of constructing the baryon wave functions reduces to the construction 
of the functions of A which transform under left and right rotations as the tensors of 
the irreducible representations T and J respectively. In the matrix representation 
that means: 
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exp(-iu, T ) ^ rj \. (A) = 4^ (A) ( 
^ exp(iu, J a ) V (rj) ] t . (A) = 4f J) (A) D$ ( e^ A » ) 



(A12) 



Replacing LHS in Eqs. (|A12| ) by RHS of Eqs. (|All|) and making subbstitution: A — > 1 
and exp(iu; a A a ) — > A one gets: 



4J J) (A) = c t/j Dg\^),4j J \A) = c tf D^(A) 

where c t j = ip^J J \l). 

In order to calculate the constants c t j let us observe that: 



(A13) 



4J J \^AR) = 47f> (A) Dg\L) D$(R) 



(A14) 



We will now make use of the following identity conecting elements of the Wigner 
matrices in representation T and its complex conjugate T: 



(A15) 



where Q denotes the charge of the state: Y, I, J 3 and minus before the index i or 
j denotes the complex conjugate state: —Y,I,—I 3 . Putting A = 1 in Eq. ([A14|) we 
we find that ip(lJR) is proportional to the product of D^ r \L^) D^(R). On the 
other hand, following Eq.( A13|) , tp(lJR) should be proportional either to D^(L^R) 
or D^\L^R). That means immediately that 



T = J. 



One can also easily convince oneself that: 



c 't.i 



So we get two equivalent expressions for the baryon wave function: 



4J\A) = (-)«Ci) D W {A t) = (_)«<*) D (J) (A ). 



(A16) 
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Throughout this derivation we have assumed that charges are integer. Finally making 
use of the fact that D^ ] {A^) = D { J\A)* we get our final expression for the baryon 



wave function: 

V(Y T T 3 ) (Y'J J 3 ) 



(A) = VdimT (- 



\Y'/2+J 3 



<Y,T,T 3 \D^ T \A)\-Y',J,-J 3 > 



(A17) 



where the normalization factor has been included. Remember that the right hyper- 
charge is constrained: Y' = — iV c /3 

The action of the colletive operators on these wave functions is straightforward: 
flavor operators T a and spin oerators Jf, (also for b = 4 ... 8) act on %j) as the SU(3) 
generators in representation T act on the state (Y, T, T 3 ) and (Y 1 , J, J 3 ) respectively. 
The action of the D functions entering the collective operators can be calculated with 
the help of the SU(3) Clebsch-Gordan coefficients: 



dim T 3 J dA D™(A)* D%f 2 {A) Dgg(A) = £ 



t\ ^2 ^3 



/ 



% T 2 V 

ji h 33 



, (A18) 



where 7 is the degeneracy index. In Tab. |XI| we list some of the proton spin up matrix 
elements of the collective operators which enter the expressions for g\. 



APPENDIX B: DERIVATION OF THE REGULARIZATION FUNCTIONS 

FOR THE AXIAL CURRENT 



Here we want to give an explicit derivation of the Q° and Q 1 contributions to the 
axial current. We emphasize the method of regularization of non-anomalous quanti- 
ties from the explicit time- ordering of the collective operators within the proper-time 
regularization scheme |30|. Then the real part can be written as: 



Re Af(x) 



1 6 

25s?(x 



■Sp(to) 



du 



u 



(j>(u) exp (—uD^D) 



(Bl) 
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where 



D = d t + H + iQ E - is a imi$ A^ X a A (B2) 



and 



rf = -dt + H- iO E - is a ilAlilh tf\ a A, (B3) 



such that 



D*D = -d 2 + H 2 + Q 2 E - i[fi E , F] - *{ft E , <9J - zs^TiTiTs^AaA, if} 

+s?747*75[A t A a A, n E ] - 2747,7 5 [s^ f A a A, (B4) 

Then one has to expand D^D around DqDq with Dq = dt + H, i.e. one expands in 
terms of s" and fi E . This is done by using the Schwinger Dyson formula: 

exp (-uD ] D) = cxp {-uDIDq) 
i 

-u J da exp (—uolDqDq) [D^D - dJd ] exp (-«(! - a)DlD ) 



o 

!. I— 

-a-P)DlDo 



W J dp J dae- uaD ° D "[D^D - DlD ]e- u P D ° Do [DW - DlD }e- u{1 - 
o o 

+ .... (B5) 
In the lowest order Q°, one obtains 



lK ' 2 5s" 



f (X%L f 

^7-^Sp ( to) / — 4>{u) / da exp (—uolDqDq) 

1 \ X ) J U J 

(uis«{l4lil5Ai\ a A, H}) exp {-u{l - a)DlD ), (B6) 



which after some simple manipulations gives eq. (^) (see Ref. ||22|| ). 

Now we want to consider the fl E corrections to the current. Let us define 
V 1 = -t[Q E ,H}-t{Q E ,d t }, V 2 = -igf{ m ys^X a A,H}, V 3 = gf m ^KA, E ] 
and V4 = —ijAjijslsfA^XaAjdt]. Consistently in the order Q E one has to consider 
combinations of V\ and V2 as well as the single sum V3 + V4. Note that it is important 
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to retain sf as time-dependent in eq. (P4j), because otherwise the two terms V3 + V4 
cancel. This can be seen using [A^Aa/L, dt] = i[£Ie, A*X a A]. After some lengthy algebra 
the operator g% defined in Eq . (|38|) can be written as: 



91 



N c t dudu' 2ujE n + 2uj'E m 
TJ 1 2^2^ (u 2 + EI)(ujV + El 



exp (i(uj — u)'){t — t )) 



poo , 

du ucf)(u) exp {-u[a(uj 2 + E*) + (1 - a)(uj 2 + )] 
< n I A c I m >< m | a 3 X b | n > T( Q^(t)D ab (t ) ). 



(B7) 



Performing now the dt integration with special care to the time-ordered product 
T( Q^(t)D ab (t ) ) one gets the relation: 



eft exp (i(u - u')(t - t )) T[Q c E {t)D ab (t )] 



PP- 



UJ — UJ 

1 



- + m5{uj — uj') 



D ab (t )ni 



pp- 



UJ — UJ 



- — m5{uj — uj') 



-pp [D ab (t ], n|] + tt5(^ - «/) {D ab (t ), 

Z — Ci/ 



^E-^afe(^o) 

(B8) 



Note that after the time-ordering the angular velocities are again assumed to be time- 
independent in order to perform the / dt integration. Last term in eq. ( |B8D vanishes 
because the 5-function makes the integral in eq. (|B7| ) odd in uj. Therefore if the indices 
of Q E and D ab are such that [D ab , Q E ] = eq. (|B~8p gives identically zero. Evaluating 
the uj,uj' integration finally gives: 

N c ifcdbD ad 



9l 



< n I A c I m >< m \ a 3 X b \ n > K Q {E n ,E m ), 



(B9) 



^ ± cc m,n 

where the regularization function is given by 
K.Q{E n , E m ) = — J du (f>{u) I da p= exp (-u[aE n + (1 - a)E m \), 







a{l — a) 



(BIO) 



which in contrast to the regularization function for the usual moment of inertia 
TZx(E n , E m ) or 1Zp(E ni E m ) is antisymmetric with respect to E m and E n . The du 
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integration can be performed analytically in the case of step like regularization func- 
tions 4>{u) = cfi{u — 1/Af) and gives 



Qt " mj " 1 Jo 2n aEl + (1 - a)i£ 



Using the formula 



da 



1 



7T 



:, < p < q 



'a(l - a) Q- a P ^Jq(p-q) 
the infinite cutoff limit of eq. (|48|) is given by (p. 219 of Ref. [|59] 

1 signE n - signE m 



E n — E n 



and was used in Ref. [17| to calculate the 1/N C corrections. Defining 

Nr. 



Q 



be 



4 



J2<n\\ c \m><m \ a 3 X b \ n > U Q {E n , E m ) 



the operator g\ can be rewritten as: 

9 a — 



D 



u3 



2iQ 



15 



D 



o3 • 



(Bll) 



(B12) 



(B13) 



(B14) 



(B15) 



h h 

Equation ( B15 ) follows directly from the real part of the Euclidean effective action 
given by eq. (BI). Therefore the corrections described above have no counterpart in 
the Wess-Zumino term which follows from the imaginary part of the Euclidean action. 
As such they vanish identically in any local mesonic theory like the Skyrme model for 
instance. 



APPENDIX C: COMPARISON WITH THE GRADIENT EXPANSION 

In order to check the results of the numerical diagonalization one should always 
consult the long wave length expansion of the coefficients appearing in the expressions 
for the observables. This technique is described at length in Ref. |45j. It also clarifies 
the question, whether the exact numerical value can be approximated by the gradient 
expanded quantities, or in other words, whether the local mesonic theory like the 
Skyrme model for instance, is a good approximation to the NJL-model. 
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1. The lowest order result from the real part of the EEA 



For the lowest order (Q°) only the quantity M3, which already exists in SU(2), 
contributes to . Its gradient expansion can be found by expanding Dq D = —d 2 + 
M 2 + iMjidiU(x) in eq. ( p6|) in terms of the gradients diU(x). The result is: 

M| rd = lJd 3 x (adm - ttAct) . (CI) 

eq. (|CH] ) can be rewritten in terms of the chiral angle 9 and for tt and a on the chiral 
circle a(r) = cos9(r) and 7r(r) = sin6 l (r): 

Mt = Y%j dr r2 i 6 ' + 2shl9 r COS9 ) • (C2) 
For the simplest case of a linear profile 9{r) = 7r(l — r/2R), it reduces to: 

Mt = (l - JL) . (C3) 

This quadratic behaviour of M| rd is explicitly checked by using a large R profile func- 
tion 9 as an input for the quark wave functions of the exact formula for M 3 . 

Another quantity from the real part of the EEA emerges due to the presence of a 
finite m s . One obtains in leading order: 

R?s d = i / cfxd^ix), (C4) 

where we have used the formula for the normalization of the kinetic term for mesons, 
which fixes the cutoff for the quadraticly divergent integral I2(M). As a total diver- 
gence the behaviour of this term is determined just by the asymptotics of the pion 
field. 

Therefore it vanishes for the linear profile discussed above or in the case of finite m 7T , 
when the profile decreases exponentially. In order to check our numerical calculations 
we have used 9 = 2arctan(— R 2 /r 2 ). Then Rf 3 = —(47r/3g)f n R 2 . As stated in the 
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text, such term is not found in the Skyrme model and it would actually give a vanishing 
contribution due to the asymptotics of the profile. In the present non-local quark model 
B,83 does have a non- vanishing contribution (compare with Tab. |VT|) . 

2. The anomalous terms from the imaginary part of the EEA 

The axial vector current gets a contribution from the imaginary part of the EEA, 
which is non- vanishing only in the SU(3) case. In a local mesonic theory it can be 
derived from the Wess-Zumino term ||34|| . Here want to show shortly how to derive this 
contribution from the non-local EEA of the present NJL model. Consider a quantity: 



l mAt{x) = -JvA(t)— St>M 



1 1 



i liia ^ b D ab s a {x). (C5) 



with D = dt + H + iVLe- Going to the operator form and using D^D one can write 

Imi°(x) = \ 5 J( x ^ P{to) jjjp D ] i^^ i ^\ b D ab -i<y i <y i <y b \ b D ab D s a (x) (C6) 

Expanding D again in terms of the gradients leads after some laborious algebra: 

Mif = J^L eQflu3 e 3ab J d z xd^ a {x)d^\x)a{x). (C7) 

In the case of the hedgehog Ansatz i.e. cr(r) = cos^(r) and 7r(r) = sin 9{r) Eq. (|C7|) 
reduces to: 

JV P W , N, 



M -— [ drrd r a(r)ir(r) 2 = — [ drr6'(r) sin 3 0(r). (C8) 

6n J2. J 6rr J 



Maa 3 

J TT 



For the linear profile 6{r) = ir(l — r/2R) we obtain a compact expression: 

Mff = -1-R. (C9) 
This linear behaviour for large size chiral fields can be seen in Fig. 7. 
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3. The f^-terms from the real part of the EEA 

In this Appendix we derive the gradient expansion for the non-local terms (Q a b)- 
The formulae below are given only for this part of the axial vector current operator 

9a- 

lr. 1 1 

9 a = ~ o Tr w / d s x dt <x,t \ \txt\ „ 7o7»75Ab \x,t > 

l J Ot + ri Ot + ri 

T[n c (t)D ab (t )] } (CIO) 

where the regularization is neglected here for simplicity. Inserting eigenstates of d t and 
H and using eq. (|B8|) we can define 



9l 



y3 y3 
A 12 , A 45 



where the X-quantities can be calculated from: 



D a3 , (Cll) 



X c = Tr / — — PP < x , A c , -— 7oTiT5 A 6 x > . C12 

Then the recipy is to multiply denominators and numerators by the hermitian conju- 
gate of the denominators and recover H 2 = —df + M 2 + iM^iid^J[x) in denominators, 
which can be expanded in terms of the gradients. Then these expressions can be 
straightforwardly simplified to the pure SU(2) quantity: 



4o7T J J w 



N C M 

N < M r dr r 2 (e> + ^ 



J d 3 x^ (ad^ - n%a) 



\ r 
N C M r „ / , _ sin2# 



12 

and the pure SU(3) quantity 



J drr 2 6' cos 26 + ^-\ (C13) 
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N C M r „ / , sin 26 



Jdrr*(6' + S -^), (C14) 



where the first line for Xf 2 is a total divergence and vanishes for chiral fields, which van- 
ish at least as 1/r 2 for r — ► oo. Assuming physical profiles, which vanish exponentially 
with the pion mass, the axial vector current operator can be written as 

sin 29 s 



g% = J d 3 r r 2 [e' 



8vr „, M M 

— f H 1 

3 Jn Ah 8/ 2 . 



D a3 . (C15) 



Note that h,l2 ~ N c , such that the last two terms in eq. (|C15| ) represent a 1/JV C 
correction. Therefore eq. ( |C15|) resembles very much the result of Dashen and Manohar 



41ji0|, which states that the 1/N C corrections to the axial current lead only to a 
renormalization of g&. Or in other words, the ratio of different coupling constants has 
no 1/N C correction. 
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FIGURES 

FIG. 1. Constituent quark mass M as a function of of m s such that the octet-decuplet 
splitting is reproduced. 

FIG. 2. The model deviations (i.e. theory-experiment) for octet and decuplet baryons as 
functions of m s for the optimal constituent quark mass M of Fig. 1 

FIG. 3. The theoretical deviations from the experimental values of the E and A particle 
for M = 419 MeV. Compared are a perturbative treatment (pertl) of the wave functions 
in zeroth oder, a linear corrections of order m s (pert2) and an exact diagonalization (YA), 
proposed by Yabu and Ando [37]. 

FIG. 4. The coefficients c\ Q and C27 of the higher representations 10 and 27 of the proton 
as a function of the constituent quark mass M. The strange current quark mass is chosen as 
m s = imMeV 

FIG. 5. The anomalous quantity M44 compared with the leading term of a gradient 
expansion, which comes exactly from the Wess-Zumino action for SU(3) pseudoscalar fields. 
This is done for a fixed linear profile in dependence of the radius R for M = 372MeV . 

FIG. 6. The g Q A , g\ and are shown for selfconsistent chiral fields in dependence of 
the constituent quark mass. The strange current quark mass is chosen as m s = \80MeV, 
according to a best fit to the hyperon spectra. 

FIG. 7. The regularization function lZQ(E n , E m ) for the time-ordered expressions (reg) 
for fixed E n and M = 400MeV in dependence of E m . This is compared to the function 
(noreg), which is obtained in the infinite cutoff limit. 
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